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Abstract 

A new 'double perturbation' theory is presented in the framework of the kinetic 
theory of quark-gluon plasma. A solvable set of equations from the 'double pertur- 
bation' is derived and are shown to be gauge-invariant. As an example, the Landau 
damping rate for the plasmon at zero momentum is calculated and discussed. 
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In the last decade, much interest was focused on the study of non-Abelian plasma from 
transport theory It is generally believed that kinetic theory can describe correctly 
the quark-gluon plasma (QGP) physics just as what can be done by the temperature field 
theory [||, and is ready to be extended to out-of-equilibrium situations Q. It has been 
demonstrated that the hard thermal loop (HTL) in temperature field theory can be obtained 
from the QGP kinetic theory but up to now there does not exist a valid scheme for 

solving the kinetic equations and at the same time keeping the non-Abelian gauge symmetry. 

A standard perturbative approach to the kinetic equations of plasma, which can keep 
the guage symmetry, is to expand the equations in the ascending power of guage coupling 
constant Another popular method in traditional study of plasma is to expand the 
equations in powers of weak field strength |[ 0. Since these two methods are in full 
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agreement with each other in the study of electromagnetic plasma ^ due to the Abelian 
nature of the dynamics (linear dynamics), they have been widely used in the literatures for 
QGP. However, they both suffer from powerlessly overcome shortcomings: The former, as 
is well known, gives only a linear leading order, i.e. an Abelian-like contribution, instead of 
truly non-Abelian ones, while the latter breaks the non- Abelian gauge covariance badly as 
was pointed out in Ref. [|[ 

If we want to solve the above-mentioned problem in the framework of kinetic theory, we 
must pay special attention to the following two aspects: the preserve of gauge invariance in 
the perturbation process and the treatment of the nonlinearity in the equations because of 
the nonlinear (non-Abelian) nature of QGP. 

In the present letter, we will propose a new scheme both satisfying the SU(3) gauge 
symmetry and having solvability for non-Abelian kinetic theory. The basic idea of the 
scheme is: After expanding the kinetic equations in guage coupling constant g, do the 
iterative calculation in powers of field strength for the purely non-Abelian counterpart. In 
the following we will call this scheme as 'double perturbation'. 

Our aim is to provide a well-defined prescription for treating the color dielectric physics. 
In this letter, we will first present the framework of the new perturbation theory and then give 
the main results of our analysis. This includes the derivation of a general set of perturbative 
kinetic equations for non-Abelian plasma and the comments of the gauge invariance and 
solvability of these equations. An application of the new formalism in the explicit calculation 
of Landau damping rate in close-to-equilibrium QGP, which has been computed from HTL 
in field theory will be given as an example. Finally, we will show some evidence that 

our results are associated with or even go beyond the HTL approximation. A more extensive 
discussion and further details on computations will be given in a longer article |p^ . 

Let us begin with the derivation of a set of equations describing the dynamics associated 
with non-Abelian fluctuations. It is sufficient to adopt the semiclassical kinetic theory of 
QGP for studying the thermal effects. The theory consists of kinetic equations and field 
equation. The kinetic equations read as |^3j 

p^^mQ±(p-^) ± fp'^5;{i^^,(a:),g±(p,x)} = 0, (1) 

P^^^G(p, x) + |p^a;{/'^,(x), G(p, x)] = 0, (2) 

where the letters with ~ represent the corresponding operators in adjoint representation of 
SU(N), in which the generator is T^. 



The background field equation is written as 



Df^F^^.^jtix), (3) 

where covariant derivative = d'^ + igA'^{x), field potential A'^ = A^t"^, and field tensor 
Ff"" = df'A'' ~ d^A^' - igfabcT^Ai^Al. Ta is the generator of SU(N). The color currents 
including external and induced ones are — j ™* + j J,"'^ ° . 

In principle, a perturbative method is often necessary to solve nonlinear equations. Here, 
for simplicity, we take the close-to-equilibrium situation of a QGP as an example to discuss 
this method. We will show later that our method is easily extended to out-of-equilibrium 
QGP. 

Following the same philosophy in traditional kinetic theory we write the distribution 
functions Q±{x,p), G{x,p) and fields as ||, [l^ , 

A; = a;, Q±=Q''^ + SQ±, G = G°^ + 5G. (4) 

The associated density fluctuations and induced fields are random quantities, satisfying 
(SQ) = 0, (SG) = and (a) = 0, where symbol ( ) denotes the random-phase average. We 
define the field tensor corresponding to as /^j,. 

In this way we obtain a set of kinetic equations for fluctuations from equations (1) and 

(2) 

P^'D^SQiip^x) ± ^p''d;{U,{x),Q±ip,x)} = 0, (5) 
P^D,SG{p, x) + G(p, x)} = 0. (6) 

In order to ensure the consistency of soft covariant derivative, in the following we will 
impose a limitation on the amplitude of the flclds: if ^ T, then ga^ ^ gT is of the same 
order as the derivative of a slowly varying quantity, id^ ^ ga^ [ p^ . 

We now employ our 'double perturbation' scheme to expand the distribution functions in 
powers of g in the first step, and then iterate repeatedly the nonlinear parts in the expanded 
equations in field quantity. These can be expressed as 

5Q = gQ(^) + g^Qi^) + . . . ^ <5G = .gG(i)+.g2G(2)+..., (7) 

q(«) = ^ q("^-^), g(") = ^ G'"'^). (8) 
Thus we obtain a series of equations 

p^^d^Q^^^ ±igY,P^K.Q^±^^] + |pn/M.,5;Q±"'^} = 0, (9) 

A 



(10) 



X 



(11) 



where n and A represent the powers of coupUng constant and field quantity, respectively. 
The induced currents associated with each order of density fluctuations are determined by 



where Nf is the number of quark flavour. The factor 2 accounts for the spin degrees of 
freedom. 

Equations (9), (10), (11) and (12) form a basis of perturbation theory of non-AbeUan 
kinetic theory. 

Some short comments follow. 

1. The non-Abelian gauge symmetry is exactly preserved in the perturbation equations 
(9) and (10) for each order of g. As a consequence Q^J!^ and G'"', like Q± and G, transform 
separately as SU(3) gauge-invariant scalars. Because the summation of non-Abelian terms 
in (9) and (10) can be done infinitely, the 'double perturbation' method we use guarantees 
the gauge-invariance of the result of a physical quantity within any high precision (A being 
an arbitrarily large number). 

2. In general, it is a difficult task to straightforwardly solve the equations (5) and (6) 
due to the the nonlinearity of non-Abelian counterparts | [l5| |. However, the summation over 
A in equations (9) and (10) imply that the iterative procedure in the powers of field quantity 
ensures the solvability of the non-Abelian counterparts. 

3. As was pointed out in Rcf. ]l6| , the perturbative expansion in field quantity only, 
which had been carried out by some authors ^ 0, have to suffer from the disadvantage 
of breaking non-Abelian gauge symmetry at each step of approximation calculation. This 
is because the results of all higher orders in field quantity contain the contributions from 
the relevant lower orders in coupling constant g for the density fiuctuations. So, in order to 
get a gauge- invariant physical result, the problem of resummation of all the contributions 
of the same order in g has to be taken into account in the past works 0, [l6| . While our 
approach automatically singles out all the contributions of the same order in g and collects 
them together. In this sense, the treatment of non-Abelian contributions in our approach 
gives a theoretical basis for the resummation technique. 

4. It can be easily verified from Eq.(12) that not only the total color current but also 




(12) 



each order current obey the covariant conservation law 



jjt.j{n)inda (13) 



This is automatically consistent with field equation (11). 

5. Our theory can be easily generalized to the out-of-equilibriuni situations through the 
decomposion of the distribution functions Q±{x,p)^ G{x,p) and fields into regular terms 
and density fiuctuations. The expression (4) is then replaced by 

^^ = (^P + <, Q±=-{Q±)+5Q±, G={G)+SG. (14) 

Such kind of division has been used in the study of QGP |4| as well as in other works [Q, ^ 
[l6| . In particular, the gauge consistency of the decomposition has been proved in Ref. JlSf . 
Our stress here is to put forward a gauge-consistent fiuctuation dynamics(a^ ^ T), while 
Litim ct al think that the regular parts describe mean field dynamics |l^ in full accordance 
with the effective soft field dynamics ((^^) ^ gT) p7| |. 

6. The 'double perturbation' approach gives an insight into the plasma with color freedom 
of degrees. 

We will get a more penetrating understanding of these remarks from the calculations 
below. 

As an example, we now apply the above formalism to calculate the gloun damping rate 
for a purely gluonic gas in the close-to-equilibrium situation. The first order equation in 
coupling constant is 

p^9^G(i) +z.g^p^[a^,G(i'^)] + <?P^/W9;g(") = (15) 

A 

^'^(a)4l^'^=J^'^'"'"(^), (16) 
Denoting the summation term by S, we have 



(27r)4 (27r)4 (2^)4 ^ P ■ {k2 + ka)? ■ ks 

x[p-a{ki),[p-a{k2),p^fJJ{k3)]] 



(17) 



It is worth noticing that all the terms in S are of same order in g. This is because 
the factor g in the numerator and the 4-dimensional wave vector fc^ ^ g) in the 



denominator appear in pairs and cancell the factor g each other. The remaining dependency 
on (7, coming from the field strength /^Jy', is of first order. 

A general formula of G'-^^(w,k) in the momentum space can be written down from 
Eq.(15). We know that the current induced by the fluctuations can be expressed as 

J^'^'"' ^=9 I 7#?,-Tr[2r,G(i)(p,a;)]}, (18) 



Following the same method in Q, we can derive the response equation of medium from 
kinetic and field equations. As an example, we do this approximately by taking A up to 2 
and consider the case of zero momentum. We get from field equation (16) together with 
equations (15) and (18), 

-u^e'^''\u,Q){a''{uj,Q)) 



2 / "1^1" «- ! 
I\ g I — -— — :-7—aL020(iiJ—iJl—UJ2)- 



X (kh (a;^ (ci;i , ki ) ) (ai^ (lj2 , ki ) ) + Kt; (a? (cJi , ki ) ) (af (lj2 , ki ) ) 
+Ktt(a?(u;i,ki)){a?(aj2,ki))) 
+Ng' J ^((a'(fci)) - {a\k^))){a\u;,0)} 

where e'--'^^ denotes the dielectric function for a mode, a — I or t represent the longitudinal 
or transverse wave, respectively, ku, Kti, Ku are coefficient functions of Wi, ki, associated 
with the definite interactions of plasmons. 

We know that damping is connected with the imaginary parts in Eq.(19). Obviously, 
the second term of the right-hand side has no imaginary part. The damping will originate 
from the physical processes described by the first and third terms. Therefore, the Landau 
damping rate can be obtained easily: 

7(-H0)-(«^+a^)^, (20) 

with 

a" = 6 J kldkidLJidLj2S{mg - cji - LU2)[Kupipi + KtiptPi + KuPtPt) 

a" = 247r / kfdkiduiii- ; — ^pi{ki) ; ^(1 - ^ „ )Pt{ki)). 

J rUgkiLU^ TTigkiUJi 

where Ku = -^f^i^iuKti = ~^^Kti,Ktt = ZJ^'^**' Pi ^^'^ represent the longitudinal 
and transverse spectral densities, respectively. 

The forms of a*' and a'' show that two typical physical processes contribute to the Landau 
damping. With the constraint of the 6 function, The first one describes the self-coupling 



interaction of a mode, which has a good formal correspondence to the result of HTL in 
field theory |jri|. The other one represents the coUisional interaction of two plasmons or 
quasiparticles, which describes the long-range interaction in QGP. The numerical result of 
a'^ has been obtained to be 5.973, while the magnitude of a* is estimated to be of well- 
matched order with a^. 

In conclusion, we have proposed a 'double perturbation' approach and derived a series 
of perturbation equations for the non-Abelian kinetic theory. Our approach gives a new 
and deeper perspective to the perturbation theory and thereby provides a progress in the 
methodological problem in the study of non-Abelian kinetic theory. 

The study of the dielectric physics in QGP will be advanced by applying this new theory. 
This includes to get a non-Abelian gauge-consistent dielectric tensor and discuss the response 
of a QGP to external sources and so on. We would emphatically point out that 

since our approach has led A-point functions (or correlators) into first-order fluctuations in 
coupling constant g, the physics beyond an equilibrium state will play more vital role in a 
QGP than in an electromagnetic plasma 

As an example, we have studied in the last part of the present letter the Landau damping 
in the zero momentum case using this approach. The result shows that both the self-coupling 
interaction of a mode and the interaction of two plasmons contribute to the damping rate. 
The physical mechanism for two-plasmon interaction can clearly be expressed as a two-body 
'collision' from long-range interaction of two quasiparticles. We need to make it clear further 
as for the self-coupling interaction, which has showed, in form, the similar 2-2 scattering 
and 2-3 scattering processes contained in the self-energy of HTL . 

We also believe that this approach will provide more and better results to uncover the 
deeper link between kinetic theory and field theory searched by many investigators 1^, |l^ . 
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